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ABSTRACT 



Context. The origin and evolution of the intracluster medium (ICM) are still not fully understood. A better understand- 
ing is not only interesting in its own right, but it is also important for modeling hierarchical structure growth and for 
using cluster surveys to determine cosmological parameters. 

Aims. To determine if there exists any evidence for evolution in the temperature or iron abundance gradients between 
z — 0.14 and z m 0.89, therefore elucidating the origin of energy and metal input to the ICM. 

Methods. By using a sample of 35 observations of 31 clusters of galaxies found in the archival data of Chandra and 
XMM-Newton with redshifts between 0. 14 and 0.89, we derived the temperature and iron abundance radial profiles. To 
compare clusters with similar properties, the data were divided into comparable subsets. 

Results. There is no substantial evidence to suggest that the iron abundance radial profiles in galaxy clusters evolve 
with redshift in any of the chosen subsets. Temperature radial profiles also do not appear to be changing with redshift 
once selection effects are taken into account. 

Conclusions. The lack of evolution in the iron profiles is consistent with scenarios where the galaxies in clusters are 
stripped of their gas at higher redshifts. The temperature and iron abundance profiles also suggest that the primary 
source of heating in high redshift clusters is the gravitational infall of mass. These findings further emphasize the im- 
portance of modeling the local environment of clusters in cosmological studies and have important implications for 
studies that go to larger redshifts. 

Key words. X-rays: galaxies: clusters 



1. Introduction 

Many dynamic processes occur in the high temperature, 
iron rich intracluster medium (ICM) of rich clusters of 
galaxies. Radiative cooling, supernova driven winds, and 
ram pressure stripping are just a few of the processes by 
which thermal energy and iron may be injected into the 
ICM. The theor y of these processes is discussed in detail 
elsew here (e.g.. lSarazinlll988clRosati et alj|2002t lArnaudl 
I2005I) . In order to help resolve which processes most sig- 
nificantly contribute to the heating and iron enrichment 
within clusters, the location and evolution with redshift 
of both iron abundance and thermal energy can be used to 
provide significant clues. 

Both Chandra and XMM-Newton have the angular 
resolution to make two-dimensional map s for nearby 
cluste rs such as Coma and Perseus (e.g., iFabian et alj 
2006), but the statistical significance of the archival data 
is not generally sufficient for these types of measurements 
in distant clusters. 

Many have already investigated how overall temper- 
atures o r iron abundances m ig ht vary with redshift, in - 
cluding iBalestra et alj d2007l). iBonamente et ail d2006l) . 
iTozzi et alj d2003l) . and [Sanderson et al.l d2006l) . So far 



very little evidence for the evolution of the ICM in 
galaxy clusters up to z — 1.0 has been discovered 
dMatsumoto et alj[2000l) . but even subtle hints of changes 
with redshift could have a profound impact on the general 
understanding of the processes involved in the ICM. This 
current study spatially resolves the temperature and metal 
(iron) abundance to higher redshifts than previous works, 
as well as determines whether there exists evidence for 
changes in the temperature or iron abundance profile with 
redshift. 



For higher redshift clusters, studies of the evolution 
of the radial temperature and metalli city profiles has been 
reported out to a redshift ~ 0.3 dLeccardi & Molendl 
2008b®, and here the work is extended out to ~ 0.9. 
The goal is to search for evolutionary changes in the ra- 
dial profiles. The results can be used to determine the 
relative significance of processes such as ram pressure 
stripping, Active Galactic Nucleus (AGN) activity, and 
galactic winds in clusters of galaxies. The end results 
can then be compared directly to simulations such as 
iDomainko et alj d2006l) . Throughout this paper the con- 
cordance model of cosmology with Hq =71 km s Mpc~' 
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Ao = 0.73 and £2 m o=0.27 is assumed. A ll calculations of 
cosmic distances and times are based on lWrighll d2006l) . 



net exposure times after all processing are listed in Table 

m 



2. The Data 

All data in this study were taken from the Chandra and 
XMM-Newton public archives, and are listed in Tables[T]& 
[2] Most o f these observatio n s were chosen after being in- 
cluded in iBonamente et al.l d2006l) . iTozzi et al.1 (120031) or 



iLeccardi & Molendil d2008b i All clusters listed here that 
are not found in these studies were found by searching 
the Chandra & XMM-Newton archives for galaxy cluster 
observations. The resulting sample was chosen to cover 
a large range of redshift with reasonable uniformity. This 
allowed for the sample to be conveniently divided into 
subsets which will be compared to on e another, build- 
ing o n previous work at lower redshifts dDe Grandi et al.l 
2004). The redshifts and coordinates of these clusters 
were provided by the NAS A/IPAC Extragalactic Database 
(NED), and the hydrogen column densities are given 
by the n# calculator provided on the HEASARC web- 
site. Investigations into the element abundances of galaxy 
clusters suggest that most of the elemental emission is 
dominated by iron, and that many other elements have 
abundance levels that do not change much with redshift 
dBaumgartner et all 120051) . This allowed for the use of 
iron abundance and metal abundance interchangeably, 
and although the model used in this work fit the abun- 
dance of all metals, the results of these fits will hereafter 
be called the iron abundance measurements. 



2.1. Chandra Data Processing 

All Chandra observations used herein were observed with 
the ACIS camera. All processing of Chandra data was 
carried out using CIAO 3.3 and all the contained pack- 
ages including Sherpa and ChlPS and fol lowed the gen- 
eral procedure found in lTozzi et al.l d2003l) . For each data 
set the standard issue EVENTS LEVEL 1 file was reset to 
remove all corrections already implemented on the data 
set. The process acisrunbadpix was run on the reset file 
to detect all bad pixels. Next, the standard Level 1 events 
processing, acisprocessevents, was run on the reset file to 
create a new EVENTS LEVEL 1 file. Detection of very 
faint pixels was done whenever the data itself was taken in 
VFAINT mode. The new Level 1 file was then filtered to 
only include the standard event grades 0,2,3,4, and 6 and 
also filtered further with the pipeline good time intervals 
to create a new EVENTS LEVEL 2 file. The new Level 2 
file was processed further using the destreak routine and 
then using the analyzeltcrv routine as part of the ChlPS 
package to determine in greater detail the good time inter- 
vals. The bin time used for all data sets was 200 seconds. 
The de-streaked file was then filtered again using these 
good time intervals. As a final step in processing, the im- 
age was then filtered only to include the energy range of 
0.3-10 keV. Filtering on the energy improves the signal 
to noise ratio of the image. The resulting event list was 
then used in all subsequent imaging and spectral analysis. 
There are 25 clusters observed using 27 Chandra observa- 
tions, with redshifts ranging from 0.142 to 0.890. All of 
the cluster data sets, along with their Equatorial coordi- 
nates, redshifts, galactic hydrogen column densities, and 



2.2. XMM-Newton Data Processing 

All XMM-Newton data analyzed in this project used all 
three EPIC cameras: MOS1, MOS2, and the PN cam- 
era. The software used was SAS version 7.0, the standard 
XMM-Newton processing software. All data sets were 
processed using the images script available on the XMM- 
Newton website. The images script first runs the funda- 
mental SAS algorithms cifbuild, odfingest, epchain, and 
emchain. It then processes the standard event list further 
to improve signal-to-noise. In order to minimize the sig- 
nificance of flare events, the images script allows for the 
user to input desired Good-Time-Interval (GTI) for each 
data set. In this case the GTI intervals were chosen as 
time windows where the rate was less than 35 counts per 
100 seconds for the MOS cameras, and 40 counts per 100 
seconds for the PN camera. After GTI filtering the im- 
ages script then cleans for bad pixels based on a sepa- 
rate input script, with default bad pixels provided by the 
XMM-Newton data center. The standard bad pixel table 
provided by the images script was used here. Although the 
images script was originally designed to create high qual- 
ity images from all three cameras, the resulting event lists 
are still useful for spectral studies. Unlike the case of the 
Chandra data sets, the event lists used for XMM-Newton 
analysis have not been subject to any energy filtering be- 
fore determining radial profiles. However, all spectra for 
both instruments are restricted to the energy range from 
0.5-8.0 keV. As a final restriction, the PATTERNS on the 
events were limited to single, double, triple, and quadru- 
ple events for the MOS cameras (PATTERN < 12) and 
limited to single and double events for the PN camera 
(PATTERN < 4). Eight XMM-Newton data sets were con- 
sidered in this study. They are listed along with Equatorial 
coordinates, redshifts, hydrogen column density, and net 
exposure times after all processing in Table [2] 



2.3. Scale Lengths: The Core and Counts Radius 

In order to average clusters together in a meaningful way, 
a scale length was needed. Two scale lengths were cal- 
culated for each cluster and used to divide the cluster 
into four annular regions: the first scale length was cal- 
culated by fitting the net intensity to a one-dimensional 
beta fit and the second by dividing the cluster into four 
regions with roughly equal numbers of net counts. The 
outer radii in these regions are labeled n, r%, r^, and 
r\ in Table [3] These radial profiles were determined us- 
ing either the surface brightness (core) or the net counts 
(counts), and the background region for these calculations 
was always an annular region at least 5" wide. This an- 
nular region was usually sufficiently small to minimize 
contamination from line-of-sight point sources while also 
containing enough counts (^ 100-1000) to be statistically 
significant(~ 3 - 10%). In some cases, the background re- 
gion needed to be larger to reach this level of significance, 
but a region with a significant number of counts ( > 100) 
was always used. The profiles are not strongly sensitive 
to the number of counts used in the background regions, 
which were chosen specifically to be radially symmetric. 
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Table 1. Basic information about Chandra data sets used in this study. 



Clusters analyzed from Chandra 



Cluster Name 


Obs # 


RA 


DEC 


7 


777777 7\ 

n H (10 i2 cm 2 ) 


Exposure Time (s) 


Abell 1413 


1661 


11 55 18.20 


+23 24 28.80 


0.142 


0.0219 


9749 


Abell 2204 


499 


16 32 47.00 


+05 34 33.00 


0.152 


0.0567 


11250 


Abell 665 


3586 


08 30 45.20 


+65 52 55.00 


0.182 


0.0431 


28600 


RX J0439.0+0520 


527 


04 39 02.20 


+05 20 43.00 


0.208 


0.1070 


10830 


Abell 773 


533 


09 17 59.40 


+51 42 23.00 


0.217 


0.0126 


10500 


Abell 697 


4217 


08 42 53.30 


+36 20 12.00 


0.282 


0.0341 


18880 


Abell 611 


3194 


08 00 56.90 


+36 03 26.00 


0.288 


0.0499 


35990 


MS 1008.1-1224 


926 


10 10 32.33 


+12 39 32.18 


0.306 


0.0726 


42960 


MS 2137.3-2353 


928 


21 40 12.70 


-23 39 27.00 


0.313 


0.0355 


31990 


Abell 1995 


906 


14 52 50.40 


+58 02 48.00 


0.319 


0.0145 


43830 


ZWCL 1358+6245 


516 


13 59 50.60 


+62 31 04.00 


0.328 


0.0193 


45420 


MACS J2228.5+2036 


3285 


22 28 34.40 


+20 36 47.00 


0.412 


0.0429 


19700 


MACS J2214.9-1359 


3259 


22 14 57.30 


-14 00 14.00 


0.483 


0.0328 


17020 


MACS J1311.0-0310 


3258 


13 11 01.70 


-03 10 41.00 


0.490 


0.0188 


15000 


MS 0015.9+1609 


520 


00 18 33.86 


+ 16 26 07.75 


0.541 


0.0407 


66940 


MACS 1423.8+2404 


4195 


14 23 47.80 


+24 04 41.40 


0.545 


0.0238 


113400 


MS 0451.6-0305 


529 


04 54 10.90 


-03 01 07.20 


0.550 


0.0500 


16370 


MS 0451.6-0305 


902 


04 54 10.90 


-03 01 07.20 


0.550 


0.0500 


41470 


MACS J2129.4-0741 


3199 


21 29 26.20 


-07 41 27.00 


0.570 


0.0484 


17690 


MS 2053.7-0449 


551 


20 56 22.40 


-04 37 43.00 


0.583 


0.0462 


43100 


MS 2053.7-0449 


1667 


20 56 22.40 


-04 37 43.00 


0.583 


0.0462 


43900 


MACS J0647.7+7015 


3196 


06 47 50.20 


+70 14 55.00 


0.584 


0.0563 


18850 


CL 1120+4318 


5771 


11 20 07.60 


+43 18 07.00 


0.600 


0.0208 


19740 


MACS J0744.8+3927 


3197 


07 44 53.00 


+39 27 26.00 


0.686 


0.0568 


19690 


MS 1137.5+6625 


536 


11 40 23.30 


+66 07 09.00 


0.782 


0.0121 


119500 


RX J1716.9+6708 


548 


17 16 52.30 


+67 08 31.20 


0.813 


0.0372 


50350 


CL J1226.9+3332 


3180 


12 26 58.20 


+33 32 48.00 


0.890 


0.0137 


29470 



Table 2. Basic information about XMM-Newton data sets used in analysis. 



Clusters Analyzed from XMM-Newton 



Cluster Name 


Obs# 


RA 


DEC 


z 


n H (10 i2 cnr 


2 ) Exposure Time (s) 


Abell 1763 


0084230901 


13 35 17.20 


+40 59 58.00 


0.223 


0.0082 


19500/19500/9197 


Abell 2390 


0111270101 


21 53 34.60 


+ 17 40 11.00 


0.228 


0.0619 


13900/13900/9445 


Abell 1835 


0147330201 


14 01 02.00 


+02 51 32.00 


0.253 


0.0204 


78120/79450/32060 


RX J0256.5+0006 


0056020301 


02 56 29.51 


+00 05 28.70 


0.360 


0.0650 


19190/19170/7914 


RX J03 18.2-0301 


0056022201 


03 18 28.76 


-03 00 46.70 


0.370 


0.0505 


20270/20340/11570 


RXJ0426. 1+1655 


0056020401 


04 26 04.20 


+ 16 55 48.50 


0.380 


0.1910 


18820/18820/8337 


RXJ1120.1+4318 


0107860201 


11 20 00.91 


+43 18 15.10 


0.600 


0.0203 


22140/22150/16240 


CL J1226. 9+3332 


0200340101 


12 26 58.00 


+33 32 54.00 


0.890 


0.0138 


77750/77780/60660 


Note: The 


exposure times 


for the XMM-Newton data 


sets are 


given in the 


following order: 



MOS1/MOS2/PN 



For each cluster, both radii have been calculated using an- 
nular regions 2. "5 wide to determine the intensity profile 
across the cluster. Both scale radii were then used to deter- 
mine regions for extracting spectra and measuring the av- 
erage temperature and iron abundance profiles. The only 
profiles and measurements shown in the text will use the 
counts radius, the reasons for which will be discussed in 
32.3.31 The measured core radius and counts regions are 
found in Table [3] 

2.3.1. The Core Radius 

The basis for using a core radius length scale is that clus- 
ter dynamics and the evolution o f the ICM could be inter- 
related (e.g., iRosati et~ai1 |2002|) . so a one-dimensional 
beta fit was performed on each cluster to determine the 
core radius scale. The centers were chosen by eye, and 
consistent with all attempts to determine the center of 
emission analytically. The error for choosing this center 



was always a factor of 10 smaller than the calculated core 
radius. The radial profile of surface brightness was calcu- 
lated using the the CIAO command dmextract, and the 
net counts annuli described in 32.31 and it was fit to the 
Sherpa betald model. The betald mode fits the data to a 
function of the form 



I{x) - Ax 



1 + 



-3/3+1/2 



(1) 



where r c is the core radius. The parameter /3 is a mea- 
sure of the ratio of kinetic energy in the galaxies mov- 
ing i n the cluster to th e thermal energy of the cluster 
(e.g., IRosati et alJl2002l) . The radial fits are described by 
the two parameters r c and /3. The offset from zero (xo) 
was always fixed at zero, and the proportionality con- 
stant normalizes the fit to the particular data relating it 
to the central density of the cluster. Radial profiles are 
listed in Table [3] with the first eight observations being 
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from XMM-Ne wton . More elaborate f itting profiles like 
a double-/? fit (iBonamente et al.ll2006l) could have been 
used in lieu of a single-/? fit, but a single-/? fit is usually all 
that statistics require (e.g., lEttori et aT1l2004al) . For clus- 
ters with high surface brightness in their center, however, 
a more appropriate fit may be the more elaborate model 
used by iKotov & Vikhlininl d2005l) . 



2.3.2. The Counts Radius 

Upon examination of the one-dimensional beta fits and 
their reduced x 1 statistics, it appeared as though the 
beta model fits derived from Sherpa were not statisti- 
cally robust. The reduced x 1 f° r these fits was often 
above 2, and although othe rs have used similar results 
dSakelliou & Ponmanll2004l) . it was judged necessary to 
find a second method used as to of independently verify 
the profiles. Therefore, a second set of temperature and 
iron profiles were calculated with the same set of data. 
The second method used was to divide the cluster into 
radial regions not by a constant length scale, but instead 
by giving each region a constant number of counts. The 
furthest extent of spectral extraction was also determined 
by the net counts in a 2. "5 wide region. If the net counts 
next the 2. "5 wide region out ever dropped below 20 for 
Chandra or 100 for XMM-Newton data, then this defined 
the edge of the cluster due to their low signal-to-noise 
ratio of the regions beyond this point. The four calcu- 
lated counts regions for each cluster are listed in Table 
[3] Dividing the cluster in this manner also allows for all 
of the spectra for one cluster to be comparably signifi- 
cant instead of having a wide disparity in spectral quality. 
The spectra were extracted in exactly the same manner as 
the core radius spectra, which will be described in detail 
in 92.51 This process of sub-dividing the cluster will be 
called the counts profile or counts radius, hereafter. They 
will also be listed as r\, r%, r^, and r^ as an abbreviation in 
the tables, particularly in Table [3] 

2.3.3. Differences Between the Two Profiles 

The radial profiles in each cluster are compared in two 
ways in Tables |3] & H Table [3] shows the extents of each 
individual counts region as well as the calculated core ra- 
dius and virial radius 7*200 m arcseconds. Table|4]shows the 
average core and counts regions in terms of kpc. It was 
found that there were insignificant differences between 
the two radii on average, except for the very outer region. 
The outermost counts regions are usually much wider 
than 2r c (see the averages given in Table|4]i. Therefore, the 
average temperature and iron abundance measurements 
should be strongly correlated between the core and counts 
radius, with differences only being expected in the outer- 
most regions. Also included in Table [3] is the calculated 
virial radiu s in arc seconds. Th e virial radius is calculated 
as found in lJones et al.l d2003l) . divided by 1 .4 to take into 
account the di fference between c osmologies (Ho = 50 km 
s" 1 Mpc~' in iJones et ail d2003l) while H = 71 km s" 1 
Mpc~' here). The formula itself is given as 



and the temperature used was the overall temperature 
measured by the procedure described in § 12.51 This al- 
lows for a direct comparison with other work that uses the 
virial radius. Typically a counts region bin corresponds 
to about 0.1 r v i r i a i, thus the profiles discussed here extend 
out on average to about 0.4 r v i,; a i. Although the core and 
counts radius were both used as scale radii for temper- 
ature and iron abundance profiles, the measurements for 
the core radius profile have been omitted from this text. 
The core radius results have been omitted because of their 
weaker statistics as well as the presence of two very small 
core radius measurements, noted in Table [4] The one- 
dimensional yS-fit also suffers from bias in the presence 
of cool core clusters with high densities in the center. The 
radial profiles produced by the two different methods are 
very strongly correlated, and therefore the calculated re- 
sults were not sensitive to this choice. Using the counts 
radius allows for more robust statistics and systematics. 



Table 4. Average outer extent of the core and counts ra- 
dius regions in kpc. 



Outer extent of region in kpc 



Radius 


Region #1 


Region #2 


Region #3 


Region #4 


Core 


87 


173 


260 


347 


Counts 


84 


177 


314 


607 



r 200 = 2.779 x 



(— f 

UOkeV/ 



x(l +zr 1 ' 5 /V 1 1 Mpc 



(2) 



2.4. Excluded Data 

Each cluster had to satisfy two general conditions before 
included in this study. Each cluster first needed to have a 
sufficient number of counts to measure a temperature that 
was significantly (> 2<x) above zero. All of the clusters 
in this sample have more than an average of 1,000 counts 
per detector (3,000 for XMM-Newton observations, 1,000 
for Chandra observations) within two core radii, and the 
minimum luminosity of all clusters from zero to two core 
radii is 1 x 10 44 erg s . The second condition was mor- 
phological in nature, as it was necessary to have a rea- 
sonable /?-fit to the data. Any observations that showed 
obvious visual signs of recent large merger activity were 
immediately excluded. Beyond that, several sets were ex- 
cluded based on the values derived from the /?-fit. In these 
cases the best fit core radius was larger than the furthest 
extent of the radial profile. These sets often had /? values 
that were unrealistically large as well, usually well above 
P = 2. Since the fit was unreliable in both parameters, 
these clusters were not considered in this sample. Table|6] 
lists all the observations investigated, but not included, in 
this study and why they were not considered. 

2.4.1. The Unusual Case of RX J1 347.5-1 145 

The very luminous cluster RX J1347.5-1 145 was also in- 
vestiga ted using two obs ervations: Chandra observation 
#3592 dAllen et alj|2002l) and XMM-N ewton observation 
#01 12960101 dGitti & Schindlerll2005l) . Although it was 
sufficiently luminous and symmetric to be included here; 
the temperatures measured between these two observa- 
tions did not ag ree within 3cr. After this work began, 
lOta et al.l d2008l) found a hot bubble in the southeast re- 
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Table 3. Radial profile parameters for all clusters. The XMM-Newton clusters are listed first. 



Radial Profile Parameters* 



Cluster Name 


RA (Center) 


DEC (Center) 


P 


Core Radius 


'i 


— — — 


r i 


r A 


^"200 


Abell 1763 


13:35:18 


+40 


59 


59 


0.95 ± 0.03 


77.6™ 


32.5 


22.5 


30 


37.5 


717.0 


Abell 2390 


21:53:37 


+ 17 


41 


40 


0.572 ± 0.004 


23.4 + 0.4 


22.5 


22.5 


30 


70 


729.1 


Abell 1835 


14:01:02 


+02 


52 


40 


0.635 ± 0.002 


18.88 + 0.1 


12.5 


12.5 


17.5 


57.5 


549.8 


RX J0256.5+0006 


02:56:34 


+00 


06 


25 




49 7+20.1 
42.3 + | y 3 5 


22.5 


15 


17.5 


20 


249.3 


RX J03 18.2-0301 


03:18:34 


-03 


02 


: 59 




20 


15 


17.5 


22.5 


253.8 


RXJ0426. 1 + 1655 


04:26:08 


+ 16 


55 


14 


0.61 ± 0.02 


15.9 + 1.3 


12.5 


7.5 


10 


25 


234.9 


RXJ1120.1+4318 


11:20:08 


+43 


18 


: 05 


0.93 ± 0.07 


40.2 ± 2.9 


17.5 


10 


15 


20 


149.2 


CLJ1226.9+3332 


12:26:58 


+33 


32 


: 45 


0.85 ± 0.03 


23.6 ± 0.9 


12.5 


7.5 


10 


22.5 


146.8 


Abell 1413 


11:55:18 


+23 


24 


13 


0.85 ± 0.05 


45.3 + 3.3 


25 


17.5 


25 


57.5 


764.4 


Abell 2204+ 


16:32:47 


+05 


34 


: 27 


0.551 +0.003 


7.5 + 0.2 


10 


12.5 


22.5 


80 


558.4 


Abell 665 


08:30:59 


+65 


50 


: 12 


0.63 ± 0.02 


36.6 ± 2.2 


27.5 


20 


25 


50 


893.3 


RX J0439.0+0520 


04:39:02 


+05 


20 


: 46 


0.58 ± 0.02 


10.3 + 1.4 


10 


12.5 


20 


50 


369.4 


Abell 773 


09:17:53 


+51 


43 


47 


1.51 ±0.25 


93.2+11.3 


27.5 


15 


20 


45 


743.7 


Abell 611 


08:00:56 


+36 


03 


: 27 


0.85 ± 0.04 


28.8 + 1.4 


15 


12.5 


17.5 


55 


360.5 


Abell 697 


08:42:53 


+36 


20 


12 


0.98 ± 0.06 


67.6 + 3.8 


30 


20 


25 


67.5 


448.1 


MS 1008.1-1224 


10:10:33 


-12 


39 


: 59 


0.69 ± 0.04 


38.9 + 3.6 


25 


20 


25 


55 


328.2 


MS 2137.3-2353 


21:40:16 


-23 


39 


: 43 


0.70 + 0.01 


9.8 + 0.3 


7.5 


7.5 


15 


70 


270.5 


Abell 1995 


14:52:50 


+58 


02 


48 


1.43 + 0.15 


69.1+5.5 


25 


15 


20 


50 


382.0 


ZWCL 1358+6245 


13:59:51 


+62 


31 


: 04 


0.57 + 0.01 


13.4 + 0.8 


12.5 


15 


20 


42.5 


289.5 


MACS J2228.5+2036 


22:28:34 


+20 


37 


: 23 


1.47 + 0.29 


50.5+7.1 


15 


10 


12.5 


27.5 


278.5 


MACS J2214.9-1359 


22:14:58 


-14 


00 


15 


1.00 + 0.13 


36.6 + 5.1 


15 


12.5 


17.5 


35 


244.2 


MACS J131 1.0-0310 


13:11:01 


-03 


10 


: 31 


1.13 + 0.20 


23.7 + 4.4 


10 


10 


10 


37.5 


201.7 


MS 0015.9+1609 


00:18:33 


+ 16 


26 


07 


1.08 + 0.09 


54.2 + 4.0 


22.5 


15 


20 


57.5 


221.0 


MACS 1423.8+2404 f 


14:23:49 


+24 


04 


: 35 


0.62 + 0.01 


6.6 + 0.3 


5 


7.5 


17.5 


47.5 


162.5 


MS 0451.6-0305(#529) 


04:54:12 


-03 


00 


. 58 


1.18 + 0.19 


51.7 + 7.5 


17.5 


12.5 


15 


30 


208.2 


MS 0451.6-0305(#902) 


04:54:11 


-03 


00 


: 48 


0.83 ± 0.03 


36.1 ± 1.5 


20 


15 


20 


60 


224.6 


MACS J2129.4-0741 


21:29:27 


-07 


41 


26 


0.66 ± 0.04 


18.6 + 2.2 


12.5 


10 


15 


37.5 


185.2 


MS 2053.7-0449(both) 


20:56:22 


-04 


37 


: 47 


0.66 ± 0.05 


18.2 + 2.1 


12.5 


10 


12.5 


32.5 


208.6 


MACS J0647. 7+7015 


06:47:50 


+70 


14 


: 55 


0.91+0.08 


29.8 + 3.1 


15 


10 


15 


37.5 


233.8 


RXJ1120.1+4318 


11:20:06 


+43 


18 


06 


1.17 + 0.31 


42.0 ± 9.7 


15 


10 


12.5 


22.5 


141.3 


MACS J0744.8+3927 f 


07:44:52 


+39 


27 


: 29 


0.56 ± 0.03 


8.0+1.4 


7.5 


7.5 


12.5 


35 


143.0 


MS 1137.5+6625 


11:40:23 


+66 


08 


20 


0.86 ± 0.07 


18.8 + 2.0 


10 


7.5 


10 


27.5 


133.1 


RXJ1716.9+6708 


17:16:49 


+67 


08 


: 27 


0.89 + 0.19 


24.4 + 5.8 


12.5 


10 


12.5 


37.5 


124.4 


CLJ1226.9+3332 


12:26:58 


+33 


32 


48 


1.29 + 0.30 


27.6 + 5.5 


10 


7.5 


10 


27.5 


141.8 



* All measurements of radial regions are given in arcseconds. t Denotes the two clusters with extremely small core radii. See ? f2~3~3l 



Table 5. Observations investigated and not included in this study, along with the reasons they were excluded. 



Clusters Not Analyzed 



Cluster Name 


Observatory 


Observation ID # 


Reason for Exclusion 


Abell 2163 


Chandra 


1653 


Bad/3-fit 


Abell 2218 


Chandra 


1454 


Bad/3-fit 


RX J122+4918 


Chandra 


1661 


Bad/3-fit 


MACS Jl 149.5+2223 


Chandra 


1656 


Bad/3-fit 


Abell 68 


Chandra 


3250 


Bad/J-fit 


RCS J0439-2904 


Chandra 


3577 


Not enough counts 


Abell 370 


Chandra 


515 


Bad/3-fit 


RX J1347.5-1 145 


Chandra 


3592 


See Section|2.4.1| 


RX J1200. 8-0328 


XMM-Newton 


0056020701 


Not enough counts 


WARP J0152.7-1357 


XMM-Newton 


0109540101 


Merger 


RX J1334.3+5030 


XMM-Newton 


0111160101 


Merger 


Sharc-2 


XMM-Newton 


0111160201 


Not enough counts 


RX J1354.3-0222 


XMM-Newton 


0112250101 


Not enough counts 


RX J1347.5-1 145 


XMM-Newton 


0112960101 


See Section|2.4.1| 


MS 1208.71+3928 


XMM-Newton 


0112190201 


Not enough counts 



gion of this cluster. Thus, it was judged that this cluster is 
not relaxed enough to be considered in this study and was 
also excluded. 



2.5. The Spectra 

Spectra for each data set were extracted with routines in 
CIAO for Chandra or S AS for XMM-Newton . If the data 
set was from Chandra then the spectra were extracted 
using the specextract routine. If the data set was from 
XMM-Newton , then the OGIP Spectral Products routine 
found in the Graphical Interface of SAS after running 
XMMSelect was used. Spectra were taken of all four an- 



nular regions in both scale lengths. When determining the 
average cluster temperature and iron abundance a circu- 
lar region from - 2r c region was used. Background re- 
gions were always circles with radii comparable or larger 
than two core radii outside of the detectable emission and 
always on the same chip as the cluster image. The back- 
ground regions remained consistent within each data set. 
These background regions are specifically chosen to con- 
tain as many counts as possible without including sources. 
Since radial symmetry is not a major concern with the 
spectral background, while the total number of counts 
is, this region is not the same background region used 
to determine the radial profiles. The software used to do 
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the spectral fitting was XSPEC 12.4.0, and the MEKAL 
model in XSPEC was always used along with the TBABS 
galactic absorption model. Since the spectra often had low 
numbers of counts, the modified Cash statistic was always 
used to determine the best-fit temperature and iron abun- 
dance simultaneously. The modified Cash statistic is ideal 
for fitting spectra with a l ow number of counts in each 
bin dNousek & Shudll989l) and also allows for the use of 
a local background spectrum instead of fitting the back- 
ground to a model. The only remaining free parameter in 
the fit was the normalization. The other parameters of the 
model were fixed: the redshifts were frozen at the values 
given by the NASA/IPAC Extragalactic Database (NED) 
while the hydrogen column densities were frozen to the 
values given by the n# calculator on the HEAS ARC web- 
site. The solar abundance values used were those from 
lAnders & Grevessd dl989l) . For the XMM-Newton spec- 
tra, all three instruments were fit simultaneously to the 
same spectral model to ensure the best use of the avail- 
able statistics. The measured temperature and iron abun- 
dance profiles for each cluster are found in Table [6] For 
the four clusters in this study with multiple observations, 
the temperature and iron abundance were determined by 
fitting all of the observations simultaneously. 



3. Selection Effects 

In any population study, it is important to take into ac- 
count selection effects within the sample itself. Two se- 
lection effects are important in this study: overall temper- 
ature and cooling time. These selection effects combined 
with the separation by redshift lead to eight different sub- 
sets of the sample which were then compared to one an- 
other. In Table [7] 8 subsets defined by their temperatures 
and cooling time are listed along with the number of clus- 
ters in each subset. 



3.1. Temperature and Luminosity 

The average temperature of this sample of clusters fol- 
lows a trend that needs to be taken into consideration. 
There is an inherent expectation that clusters at higher 
redshift should exhibit hi gher temperatur e s due to the 
L x - T relationship (e g., lOta et al.l l2006t iPacaud et all 
I2007t lEttori et al.l 120021, and references therein). This is 
because more distant (higher redshift) clusters need to be 
more intrinsically luminous to be detected than low lumi- 
nosity clusters. The X-ray bolometric luminosi ty is pro- 
portio nal to temperature as L x ,boi K T 2 5 (e.g., lOta et al.l 
[2006b . Thus on average, the observed clusters at higher 
redshift should tend to higher average temperatures. It 
was necessary to take into account this selection effect 
when searching for variations with redshift. Based on the 
available sample, the boundary between high and low 
temperature clusters was set to 6.8 keV. Setting the bound- 
ary to this temperature allows for comparably sized sub- 
sets for high and low temperature clusters, even after the 
two other selection effects (cooling time and redshift) are 
taken into account. The average radial profiles presented 
below are not strongly dependent on small changes 0.2 
keV) to this boundary. 



3.2. Central Cooling 

An important mechanism in galaxy cluster evolution is 
the process of radiative cooling. T he cooling time for a 
cluster can be well approximated by lSarazirJ (Il988h as 



f coo i = 8.5xl0 10 yr 



10 8 K 



1/2 



(3) 



where n p is the particle density in units of cm -3 and T g 
is the temperature of the gas in Kelvin. By relating the 
measured temperature and luminosity to the density and 
performing the calculation, the cooling times for each of 
these clusters can be calculated. Because high densities 
reduce the radiative cooling time, radiative cooling is po- 
tentially significant only in the central regions of a cluster. 
Therefore, all cooling calculations used the temperature 
and luminosity measured for the innermost regions. The 
cooling times were calculated assuming that the density 
is constant within the region of zero to one-half of a core 
radius and that all of these clusters originally formed at 
a redshift of z — 2. For the calculation of the density, the 
bolometric luminosity (from 0-100 keV) and the volume 
of the - 0.5r c region were used to determine an emissiv- 
ity for the region, which can then be related to the density 
bv lSarazinl (119881) as 



e s = (1.435 x 10~ 27 ) T l g /2 n\ erg cm" 3 s" 



(4) 



This assumes that the X-ray emission is all due to thermal 
bremsstrahlung and also that the density and temperature 
of the ICM is constant within 0.5r c . A single beta fit was 
used instead of a more elaborate fit, and de-projection of 
hydrost atic equilibrium models w as not done, as was car- 
ried by iBonamente et al.l d2006l) . The cooling time only 
goes as the square root of the temperature, and hence it 
can be seen that any contamination of the temperature de- 
termination by the outer portions of the cluster (in cooling 
clusters only) would be a weak effect. In comparison wit h 
the cooling times calculated in IBonamente et all d2006l) . 
there is only one discrepancy between their results and 
those derived here. The cluster MACS J2129.4-0741 has 
a cooling time slightly smaller than its age (4.64 Gyr as 
compared to 4.82 Gyr). Since this difference is small com- 
pared to the c ooling time, and i t is not listed as a cooling 
core cluster in IBonamente et all d2006l) . it is treated in this 
study as a non-cooling core cluster as well. 

The results of these calculations are listed in Table [8] 
The presence of cooling core clusters leads directly to an- 
other important process: central Active Galactic Nucleus 
(AGN) activity. For nearby clusters with cooling cores, it 
has been shown that they do no t cool as quick ly as the the- 
oretical cooling times suggest (Fabian 1994). AGN activ- 
ity could influence the temperature as far out as =; 50-100 
kpc, which is, on average, approximately 0.5r c or in the 
first counts profile bin. The majority of cool core clusters 
exhibit the radio e mission or bubbles us ually associated 
with AGN activity dDunn & Fabianl2006l) .The findings of 
iDunn & Fabianl (120061) also show that only a small frac- 
tion of cool core clusters show no evidence of AGN ac- 
tivity as either bubbles or radio sources. Therefore, it is 
expected that central AGN activity is generally (if not al- 
ways) present in the cooling core clusters in this sample 
as well. For completeness, all clusters were searched for 
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possible counterpart radio sources within one arc-minute 
of the cluster position. The presence of sources was done 
using the NVSS catalog (ICondon et all 1 19981) . and the 
sources were confirmed as being associated with the clus- 
ter using NED. The clusters with likely radio counterparts 
are denoted with a dagger in Table[8] The presence of cen- 
tral AGN activity will be assumed in discussing the results 
for cool core clusters. 



4. Results 

After dividing the sample into eight comparable subsets 
based on temperature, cooling time, and redshift, relevant 
average profiles were calculated and compared. These di- 
visions had to be made so as to make valid comparisons 
between profiles derived from low and high redshifts. 
Since there is currently little evidence for evolution in 
galaxy clusters even at high redshifts dMatsumoto et"ai] 
l2000tlMushotzkv & Loewensteinll 19971) . the null hypoth- 
esis for statistical tests will be that cluster abundance and 
temperature profiles d o not cha nge with redshift. See also 
iLeccardi & Molendil d2008albl) for z = 0.1-0.3 studies. 
Only if the statistics conclusively suggest rejecting this 
hypothesis will evolution be considered. 



4. 1 . Zero iron abundances and N/A values 

As seen in Table 6, there are several instances where the 
best fit iron abundance was zero. As a compromise be- 
tween ignoring these points altogether and producing an 
un-weighted average, the derived uncertainties to these 
points were used as statistical weights (one over the er- 
ror squared). Since these zero values were most likely 
due to an insufficient number of detected photons rather 
than a true lack of iron in the ICM, this averaging was 
judged to be a valid approach. As a cross check, these 
points were excluded from the averaging as well, and the 
results were all within 1 cr of each other in both weighted 
and unweighted averages (see 34.21 ). Also, there was one 
instance where the measurement failed to converge for 
both the temperature and iron abundance simultaneously, 
and no measurement errors could be calculated. This mea- 
surement is listed with a "N/A" in Table|6] 

4.2. Averaging Procedures 

In order to make the most robust study of the data, two av- 
eraging procedures(weighted and unweighted) were done 
on each subset listed in the following sections. The pri- 
mary averaging process used was a weighted average that 
took into account the measurement error on each obser- 
vation. For a given sample of N clusters, the weighted 
temperature/iron abundance average M is calculated us- 
ing the measurements m,- and the measurement error <x ; 
as 



M 



vTw r 



,-1/2 



2 n-2 



(5) 



As a check to ensure that no significant biases arise 
due to the weighting procedure, a standard (unweighted) 
mean and standard deviation on that mean were also cal- 
culated, both of which only depend on the measurements 



themselves and not the error on those measurements. The 
unweighted averages further demonstrate the robustness 
of these measurements, as the choice of weighted or un- 
weighted averaging does not affect the basic conclusions. 

4.3. Average Profiles 

The average temperature and iron abundance profiles are 
written out in Tables l9l&[T0l They are organized by sub- 
set and radial region, and show the unweighted values 
in parentheses next to the weighted averages. Since the 
weighted mean emphasizes measurements with the small- 
est errors and there was substantial variance in the quality 
of the spectra, the weighted mean was chosen to be the 
basis for the subsequent figures and discussion. However, 
both means exhibited the same general trends and results. 



4.4. Temperature Profiles 

The temperature profiles of Table [9] suggest that the over- 
all temperature and cooling time completely define the 
temperature profiles of clusters. The average temperature 
profiles for cooling core clusters exhibit a clear increase 
in temperature with radius, while clusters without cool 
cores are either isothermal or have temperatures that de- 
crease with radius. Similar temperature profiles for non- 
cool c ore clusters has been confirmed b y I Sanderson et af] 
J2006I) & |Leccardi & Molendil d2008bl) . 

There are also no differences between the temperature 
profiles of high and low temperature clusters of compara- 
ble cooling times with the exception of an overall offset 
due to the higher overall temperatures. Because the tem- 
perature profiles do not seem to deviate from expectations 
based on their overall temperature and cooling times, the 
discussion in f5]will focus primarily on the measured iron 
abundance profiles. 

4.5. Iron Profiles 

The iron abundance radial profiles are now taken into con- 
sideration. 



4.5.1. Low Temperature, Cool Core Clusters 

The first subsets to consider are those of low tempera- 
ture, cool core clusters. These clusters have been shown 
to hav e strong iron abundance gradients dDe Grandi et al.l 
I2004I) at low redshifts and are often associated with cen- 
tral AGN activity. Not only do the two subsets consid- 
ered here (subsets # 1 & 5) show no evidence for evo- 
lution with redshift, the iron abundance profile is almost 
identical to the profi les of still lower (z < 0.1) redshift 
clusters discussed in De Grandi et al. 1 (12004 . The results 
of iDe Grandi et al.l d2004l) have been re-projected to the 
length scale used here and corrected for the differences in 
their solar abundance model. All three profiles are shown 
in Fig.Q] From Fig.[T]there is no evidence for evolution in 
the iron abundance of cool core clusters from < z >— 0.7 
to < z >— 0. It should also be noted in Table [TO] that 
low temperature cool core clusters (subsets 1 & 5 ) do 
not have statistically significant differences in overall iron 
abundance or in average iron abundance profiles versus 
redshift. The same lack of difference is also seen when 
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12 3 4 

Radius (Counts Region) 




2 3 
Radius (Counts Region) 



Fig. 1. Iron abundance as a function of radius for clusters 
with low temperatures and short cooling times, separated 
by redshift. The dotted line corresponds to subset # 1 ( 
z < 0.4) while the dashed line corresponds to subset # 5 
(z > 0.4). The solid line is from lDe Grandi et al.ld2004l) . 




2 3 
Radius (Counts Region) 



Fig. 2. Iron abundance as a function of radius for clusters 
with high temperatures and short cooling times, separated 
by redshift. The solid line corresponds to subset # 2 (z < 
0.4). The dashed line corresponds to subset # 6 ( z > 0.4). 



comparing high temperature cool core clusters (subsets 2 
&6). 



Fig. 3. Iron abundance as a function of radius for clusters 
with low temperatures and long cooling times, separated 
by redshift. The solid line corresponds to subset # 3 (z < 
0.4 ). The dashed line corresponds to subset # 7 (z > 0.4). 




2 3 4 

Radius (Counts Region) 

Fig. 4. Iron abundance as a function of radius for high 
temperature clusters with long cooling times, separated 
by redshift. The solid line corresponds to subset # 4 ( z < 
0.4). The dashed line corresponds to subset # 8 (z > 0.4). 

4.5.3. Low Temperature, Non-Cool Core Clusters 

Subsets # 3 & 7 describe the average profiles for low 
temperature clusters without cool cores. Although both 
of these samples are small, Figure [3] shows that the two 
profiles are very consistent with one another. 



4.5.2. High Temperature, Cool Core Clusters 

There are only two clusters in the high redshift subset for 
clusters with high temperatures and cool cores (subset # 
6) which is too small of a sample to be considered repre- 
sentative of the entire population. However, it should nev- 
ertheless be noted that the average iron abundance profile 
for this small subset is entirely consistent point-for-point 
with the results of the larger low redshift sample. As can 
be seen in Figure [2] the average profile in both cases has 
a higher central abundance, with evidence for a decreas- 
ing gradient outward. Even though the average profile for 
the high redshift sample may not represent the population, 
this small sample still shows no evidence for evolution in 
the iron abundance profile with redshift. 



4.5.4. High Temperature, Non-Cool Core Clusters 

Both subsets in this case have a relatively high number of 
clusters. Again in these two subsets, the two radial pro- 
files are wholly consistent with one another, with no rea- 
son to believe that the iron abundance in these subsets is 
changing with redshift. All of the analysis suggests that 
the iron abundance is not evolving with redshift indepen- 
dent of the temperature and cooling time. 

5. Discussion 

5.1. The Data Analysis 

Before discussing the radial dependencies of the iron 
abundance and the temperature, it is worthwhile to review 
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why and how the data from the clusters were averaged to- 
gether in the manner that was done. As noted in § 2.3 and 
2.4, two different radii/scales were chosen. The first was 
the core radius derived from a single /3 model. This radius 
corresponds directly to the dynamical scale of the cluster. 
This is the natural choice for a scale length here, as it is 
based on the dynamics of the cluster. 

A second scale length was used to find four radial re- 
gions with approximately equal numbers of net counts. 
Although this scale is independent of the dynamics, it 
is useful for testing whether the profiles are sensitive to 
the scale length chosen. This "equal counts per enclosed 
area" process was chosen for its simplicity in calculation, 
and for its ability to increase the significance of mea- 
surements in the outer regions of the clusters. When the 
counts and core radii were converted into a linear pro- 
jected length scale such as Mpc, they cover very simi- 
lar ranges. The main inconsistency occurs in the outer- 
most regions, where the counts radius extends much fur- 
ther outward than the core radius. Therefore, overall con- 
sistency between the results derived by the two differ- 
ent methods was expected. Since the counts radius is so 
strongly correlated with the core radius and offers better 
statistics across each cluster, the results of the core ra- 
dius measurements have been omitted from the text even 
though they are wholly consistent with these results. 

Another possible choice of length scale is the virial 
radius of the cluster, but the easiest way of estimating 
r v iriai (or r2oo) relates d irectly to the overall temperature 
(e.g.. iJones et alj|2003l) . If the virial radius was used in 
this case, the temperature would be correlated with the 
radial regions. B ecause of the L x ,boi _ T relationship (e.g. 
lOta et al]|2006l) and clusters at higher redshifts must be 
on average more luminous to be detected with a suffi- 
cient statistical significance, the scale radius chosen needs 
to take into account this selection effect. Because the 
virial radius is related to the overall temperature, and 
therefore luminosity, it is difficult to correct for this se- 
lection effect. As a n example, there ex ist nearby clus- 
ters such as Perseus d Arnaud et al Hi 1994 and Coma (e.g. 
Watanabe et al.ll999HDe Grandi et al.ll2004tlAdami et al.l 



2006) with comparable X-ray luminosities that exhibit 
very different temperature and iron abundance profiles. 

As a further check of our analysis, our results 
we re compared with ot h er pub lished measurements such 
as iKotov & VitdTiinlnl d2005l) for CL 0016+16 (MS 
0015.9+1609), and the temperature/abundance results 
pu blished here tend to be hig her by ~ 1 - 1 ,5cr than those 
of IKotov & Vikhlininl d2005l) . The work of lBalestra etaT] 
d2007l) shows how using an updated analysis tool can 
bring out difference of this scale. Furthermore, the dif- 
ferences may be due to systematics between the Chandra 
and XMM-Newton analysis chains as well as the statis- 
tics between the two observations themselves. These dif- 
ferences in the core radius measurements are on the or- 
der of about 20% which should not strongly affect the 
results of the average pro files. The measured profiles of 
IKotov & VikhiTn in (2005]) also measure the radial pro- 
file in a diffe rent energy band t h an ou r measurement 
(0.5-2 keV in IKotov & Vikhlininl d2005l) versus 0.5-10 
keV in this wo r k). Th ree other clusters were studied by 
ISnowden et al.l d2008l) . and these results are consistent 
within 10-20% and usually within 1 standa r d devi ation. 
An important caveat is that ISnowden et al.l d2008l) used 



the MOS cameras of XMM-Newton , whereas Chandra 
was used in this study for measuring the temperature of 
Abell 1413 and Abell 2204. The systematic differences in 
response and sensitivity may be r esponsible for the dis- 
crepancies between the results of ISnowden et alj d2008l) 
and those listed here. 

The results of the counts profiles results will now be 
taken at face value, and the physical interpretation of no- 
ticeable variations will be discussed. 



5.2. Review and Comparison With Previous Work on 
the Chemical Evolution of the ICM 



Balestra et all (120071) . iMaughan et ail d2008l) . and 
Werner et al.l d2008l) all discuss evidence for the 



chemical evolution of the ICM. Older works ( e.g., 
iMushotzkv & Loewenste"ir]|l997l: iRizza et al.l 1 19981) did 
not report any statistically significant evidence for chem- 
ical evolution with redshift. The recent w ork that found 
the strongest evidence for evolution is iBalestra et al.l 
but this evolution was found only when data based on 
all the z < 1 clusters in their sample were combined 
(compare with the smaller previous set in their Fig. 3). 
Furthermore, the derived chemical evolution required the 
use of the formal statistical best fit uncertainty instead of 
the observed dispersions (e.g. their Fig. 14). In contrast, 
fo r a similar redsh ift range (below about 0.9), the results 
of lMaughan et al.l are consistent with no evolution. 

From the t heoreti c al point o f view as described in 
iKapferer etail d2007l) . ISarazinl d!988l) . and references 
therein, the chemical evolution of the ICM is a complex 
combination of effects due to cluster merging, infall of en- 
riched material, galactic ram pressure stripping, galactic 
winds, and other possibl e causes. For the relax ed clusters 
considered in this paper IKapferer et all d2007l) predicted 
no observable chemical evolution from redshift 0.15 to 
0.9. 

Our paper shows how it is important to compare like 
clusters and that there is room for further theoretical stud- 
ies to understand the apparent lack of evolution of the ra- 
dial dependence of the enriched material of the ICM, and 
at the same time induce heating via infall of metal rich 
gas clouds. 

5.3. Scenarios for heating and iron enrichment 

The data presented here are consistent with no evolution 
in radial profiles of iron abundance and temperature. The 
modeling of this should be consistent with the concept 
of hierarchical f ormation of structure in the universe (e.g. 
iGao et all2008[ and references therein). Because the iron 
abundance does not appear to be strongly correlated to the 
overall temperature, gravitational infall is the most likely 
scenario for heating clusters beyond ^6.5 keV. This in- 
fall must take place in such a manner as not to change the 
radial profile of iron abundance or temperatures. For non- 
cooling clusters this means the clusters somehow have no 
measurable temperature gradients inside about 500 kpc 
but maintain their abundance gradients while (presum- 
ably) increasing in temperature in a hierarchical growth 
model. In order to demonstrate why infall is the preferred 
major energy input, the process of heating has to be con- 
sidered in some detail. 
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An estimate of the amount of near solar abundance 
material that is added along with the energy to boost the < 
kT> fro m 4 keV to 8 keV is made now. We a s sume typical 



T dOta et al.l l2006: 



2006) and L x - M (e.g. 



Ettori et al.ll2004t ; 



Rykoffet al. 2008 



Stanek et al 



Stanek et al 



2006) relations apply to the clusters studied here. Under 
these assumptions, if the temperature increases by a fac- 
tor of 2 then the luminosity increases by at least a factor 
of 4 and the mass by a factor of 3, the amount of mass 
added to the clusters being on the order of 3 x 1O 14 M . 
Since galaxies are thought to ma ke up only about 1 /5 of 
the baryonic mass in clusters (see lLoewensteinir2006l and 
references there in), it is implausible that this mass and 
additional gravitational energy is added by normal galax- 
ies. This mass and gravitational energy must come instead 
from the infall of atypical galaxies that have gas masses 
that greatly exceed their stellar masses or gas clouds that 
never formed into galaxies. Damped Lyman a absorbers 
(DLAs ) have metallicities that are near solar at z of 1 or 
higher (Meiring et al.l2007l and references therein) which 
makes the ability to add both mass and metals via in- 
falling clouds plausible. It is beyond the scope of this 
work, though, to carry out detailed calculations of this 
infall scenario. However, an implication of the estimates 
made here is that sight lines on the outskirts of clusters 
should show the presence of DLAs at the cluster redshift. 
Furthermore, if mass is added with energy as implied by 
the L - T and L - M relationships used here, this rules 
out processes which might provide energy and metals but 
not significant amounts of additional mass, such as SN or 
AGNs. 

Instead, suppose that the L - M relationship used 
here doesn't apply to these high z clusters. In this case, 
supernovae would still not be a plausible explanation. 
This is because the energy input of approximately 10 64 
ergs would require an unreasonably high number of SNe 
when all the available en ergy is transferred to the ICM 
dConrov & Ostrikerll2008l) . For example, suppose there 
are 1,000 galaxies per cluster. This translates into 10 10 
SNe per galaxy or 10 1 SNe/year for 10 9 years. It is also 
implausible that a central AGN could provide this much 
heat (10 64 ergs), as this heat would require 10 55 erg/year 
deposited for 1 Gy, or 10 47 erg s _I minimum energy gen- 
eration assuming 100% efficiency in transfering energy to 
heat. Since the magnitude of heating seems beyond what 
central AGN activity could provide an infall scenario is a 
more likely explanation. 

5.4. Scenarios for Mixing 

Since there is no evidence for evolution in every type of 
cluster considered in this study (as seen in the figures), it 
will be assumed for the sake of discussion that clusters 
are not mixing over this range of redshifts. 

Beginning with the cool core clusters in Fig. Q] it is 
seen that they exhibit almost no signs of evolution in their 
iron abundance profiles from the 0.4-0.9 redshift bin up 
to the present day, even though there is a very clear gradi- 
ent at all redshifts. This suggests two possible scenarios: 
either enrichment and mixing both exist in such a way 
that neither is dominant, or neither process occurs. If nei- 
ther process occurs, the unchanging iron abundance pro- 
file can be explained as being due to the average clus- 
ter galaxy having lost most of its gas by z ~ 0.8. In 



this case the galaxies have no gas to stir up and mix 
the ICM. An absence of gas in cluster galaxies would 
also manifest itself in a l ow, non-evolving star formation 
rate dHomeier et al]l2005l) . A low star formation rate oc- 
curs if ram pressure stripping has removed the majority 
of the gas from most of the cluster galaxies before red- 
shift z = 0.8 within 500 kpc of the cluster center, well 
within the regions observed in this study. The suggestion 
that mo st gas is stripped from gala x ies by z — 0.4 was als o 
seen by lButcher & Oemlerl (fT98H) : lDressler et al.1 dl997l) . 

As the gradient in iron abundance has remained con- 
stant with redshift, there is also no evidence for growth 
in the gradient. Substantial changes in the abundance gra- 
dient would be expected if there was continuous activ- 
ity from the central AGN . There fore, energy input to the 
ICM beyond gravitational infall dBode et alJl2007l) must 
not have caused appreciable iron enrichment. 

For non-cool core and high temperature clusters, there 
is less statistical evidence to suggest that the iron abun- 
dance decreases with radius in the same manner as low 
temperature cool-core clusters. However, as noted above, 
it is assumed there is no evolution in the radial profiles of 
the ICM with redshift. In this case, galactic motion within 
the cluster does not result in ram pressure stripping, one 
of the key processes for iron enrichment (and possibly 
mixing) in the ICM. Without ram pressure stripping from 
galaxies to enrich or mix the ICM, the iron abundance 
(and iron abundance profile) of the ICM remains static 
over long times. 

6. Summary and Conclusions 

Average radial profiles of the temperature and iron abun- 
dance of X-ray bright clusters in the redshift range from 
about 0.14 to 0.9 were calculated. We find no evidence 
for evolution within similar sets of clusters. The total el- 
ement abundance remains constant with decreasing z and 
the gradients remain approximately the same. These ob- 
servations suggest that gravitational infall is the domi- 
nant mechanism for heating of the ICM, and also that 
ram pressure stripping does not substantially change the 
ICM out to large redshifts. If clusters experience a signif- 
icant (factor of 2 or more) amount of gravitational mass 
infall, then the iron abundance profile remains constant 
while the overall temperature increases. This infall needs 
to simultaneously keep the profile (within error of about 
20%) isothermal while also maintaining the same iron 
abundance profile. Mass and energy considerations sug- 
gest infall of metal rich gas, but how the process takes 
place needs further theoretical study. These same consid- 
erations of mass and energy suggest that although central 
AGN activity is directly related to cooling cores, the cen- 
tral AGN does not produce significant amounts of metal 
rich gas and energy to the ICM beyond its local environ- 
ment of about 50-100 kpc. 

Other theoretical challenges to consider include how 
clusters formed in the first place with both abundance 
gradients and isothermal temperature profiles. Work also 
needs to be done o n how the frequ ency and scale of 
merger activity (e.g. IWik et ai1l2008l) changes the abun- 
dance and temperature gradients of clusters. Another po- 
tentially important aspect of modeling galaxy clusters the- 
oretically appears to be an accurate modeling of the local 
environment surrounding the cluster, as the infall of mate- 
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rial initially outside of the cluster could play an important 
role in its evolution. 

From an observational point of view it will be inter- 
esting to measure temperature and abundance gradients 
for clusters at z greater than 0.9 as well as make measure- 
ments of DLAs by means of QSOs whose sight lines are 
at ^ 1 Mpc distances from the centers of rich clusters. 
These measurements may help determine more conclu- 
sively whether galaxy clusters are evolving in the ways 
described here and what kind of material might be falling 
in to cause the highest temperature clusters of galaxies. 
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Table 6. Temperature and iron abundance measurements for all clusters by counts region. 



Temperatures and Iron Abundance Measurements 



Cluster Information Temperatures in keV Iron Abundances in Zq 



Cluster Name 




Region # 1 


Region #2 


Region # 3 


Region # 4 


Region #1 


Region #2 


Region #3 


Region #4 


Abell 1413 


0.142 


6.0 + 0.4 


7.8 ± 0.7 


7.4 ± 0.7 


5.5 ± 0.4 


0.46 ±0.14 


0.41 ± 0.16 


0.49 ±0.11 


0.73 ± 0.14 


Abell 2204 


0.152 


3.9 ±0.1 


6.7 ± 0.3 


10.2 ±0.7 


12.1 ±0.8 


0.98 ± 0.09 


0.57 ± 0.09 


0.41 ±0.11 


0.25 ± 0.06 


Abell 665 


0.182 


7.7 ± 0.7 


7.9 ± 0.7 


8.5 ± 0.7 


8.1 ± 0.6 


0.32 ±0.12 


0.39 ± 0.13 


0.47 ± 0.14 


0.31 ±0.10 


RX J0439.0+0520 


0.208 


3.1 ±0.2 


4.5 ±0.3 


4.8 ±0.5 


4.2 ± 0.5 


0.50 ±0.19 


0.94 ± 0.28 


0.44 ±0.13 


+ 0.32 


Abell 773 


0.217 


7.3 ± 1.0 


9 1 +19 

-1.3 


9 6+i.s 
-i.i 


8 1+° 
-i.i 


0.43 ± 0.22 


0.27 ± 0.23 


0.68 ± 0.28 


07+ 25 

-0.07 


Abell 1763 


0.223 


7.5 ± 0.4 


8.3 ± 0.4 


7.0 ± 0.3 


6.6 ±0.3 


0.47 ± 0.08 


0.30 ±0.08 


0.30 ± 0.06 


0.26 ± 0.05 


Abell 2390 


0.228 


7.3 ±0.2 


9.4 ±0.3 


10.0 ± 0.5 


9.7 ± 0.4 


0.43 ± 0.05 


0.43 ± 0.06 


0.29 ± 0.06 


0.17 ±0.05 


Abell 1835 


0.253 


5.0 ±0.1 


6.8 ±0.1 


7.2 ±0.1 


7.4 ±0.1 


0.46 ± 0.02 


0.37 ± 0.02 


0.26 ± 0.02 


0.23 ± 0.02 


Abell 697 


0.282 


9.4 ±0.9 


9.6 ± 1.0 


9.3 ± 1.0 


9.1 ± 1.0 


0.50 ±0.16 


0.37 ±0.15 


0.54 ±0.11 


0.31 ±0.09 


Abell 611 


0.288 


6.7 ± 0.4 


7.0 ±0.5 


7.5 ±0.6 


6.5 ±0.5 


0.45 ±0.12 


0.24 ±0.10 


0.21 ±0.11 


0.34 ± 0.07 


MS 1008.1-1224 


0.306 


5.9 ±0.5 


5.8 ±0.6 


7.4 ±0.8 


6.5 ± 0.7 


0.50 ±0.14 


0.35 ±0.10 


0.37 ±0.10 


04+ 010 

u ' u -0 04 


MS 2137.3-2353 


0.313 


4.1 ±0.1 


5.1 ±0.2 


5.9 ±0.3 


5.8 ± 0.5 


0.42 ± 0.06 


0.32 ± 0.06 


0.32 ± 0.07 


0.35 ± 0.06 


Abell 1995 


0.319 


8.3 ±0.7 


9.9 ±0.8 


8.6 ±0.8 


8.8 ± 1.0 


0.33 ±0.12 




0.36 ±0.10 


0.36 ±0.11 


ZWCL 1358+6245 


0.328 


4.5 ±0.2 


8.9 ± 0.7 


8.5 ± 1.2 


9.4 ± 1.8 


0.45 ±0.10 


0.42 ±0.14 


0.27 ±0.18 


0.69 + <»° 


RX J0256.5+0006 


0.360 


5.5 ±0.4 


5.4 ±0.5 


4.4 ± 0.4 


3.7 ±0.2 


0.31 ±0.13 


0.39 + 0.14 


0.23 ±0.10 


0.38 ±0.14 


RX J0318. 2-0301 


0.370 


4.0 ±0.3 


5.5 ±0.5 


5.3 ±0.6 


6.3 ±0.8 


0.55 ±0.12 


0.27 ±0.14 


0.21 ±0.14 


0.29 ±0.17 


RXJ0426. 1 + 1655 


0.380 


4.5 ±0.4 


6.4 ± 1.0 


4.9 ± 0.7 


5.0 ±0.7 


0.47 ±0.15 




07 +016 

—0 07 


0.38 ± 0.22 


MACS J2228.5+2036 


0.412 


7.6 ± 1.3 


9.9 ±2.2 


7.9 ± 1.4 


6.9 ±0.8 


0.56 ±0.21 


0.26 ±0.19 


0.29 ±0.15 


0.52 ±0.13 


MACS J22 14.9- 1359 


0.483 


8.0 ± 1.2 


9.1 ± 1.7 


12 7 +5J 

-2.8 


8.4 ± 1.5 


0.27 ±0.21 


0.23 ± 0.22 


17+ 020 

-0 17 


13+ 018 

-0.13 


MACS J1311.0-0310 


0.490 


4.7 ±0.5 


6.8 ±0.9 


7.4 ± 1.3 


4.8 ± 0.6 


0.19 ± 0.15 


1 36+ 051 

-0.23 


09+ 20 

-0.09 


0.22 ±0.17 


MS 0015.9+1609 


0.541 


9.3 ±0.8 


9.9 ± 1.0 


10 5 +15 

Iu -0.7 


10.3 ± 1.1 


0.61 ±0.14 


0.22 ± 0.08 


0.09 ± 0.09 


0.16 ±0.09 


MACS 1423.8+2404 


0.545 


4.1 ±0.1 


5.8 ±0.2 


7.1 ±0.3 


6.9 ± 0.5 


0.66 ±0.08 


0.48 ± 0.07 


0.34 ± 0.05 


0.47 ± 0.07 


MS 0451.6-0305 


0.550 


9.9 ±0.8 


8.9 ±0.6 


9.0 ±0.7 


7.3 ±0.6 


0.64 ±0.14 


0.28 ±0.10 


0.16 ±0.10 


0.49 ±0.13 


MACS J2 129.4-0741 


0.570 


7.0± 1.1 


7.2 ± 1.2 


7.6 ± 1.2 


5.2 ± 0.7 


0.45 ± 0.24 


0.33 ± 0.20 


0.41 ± 0.23 


0.65 ±0.17 


MS 2053.7-0449 


0.583 


4.7 ± 0.7 


4.5 ± 0.6 


5.6± 1.1 


3.4 ±0.6 


0.32 ±0.24 


0.23 ± 0.20 


0.26 ± 0.23 


0.59 ±0.41 


MACS J0647.7+7015 


0.584 


14 9 +4 - 5 


11.7+IJ 


10 - 7 -2J 


9.5 ± 1.6 


0.43 ± 0.39 


0.40 ± 0.25 


0.39 ± 0.22 


0+0.18 


RX Jl 120.1+4318 


0.600 


5.4 ±0.4 


5.2 ±0.4 


5.0 ±0.4 


4.9 ± 0.5 


0.54 ±0.10 


0.51 ±0.18 


0.26 ±0.12 


+ 0.16 


MACS J0744.8+3927 


0.686 


5.3 ±0.6 


8.7 ± 1.7 


9.0 ± 1.6 


5.9 ±0.7 


97+ - 40 


0.18 ±0.23 


0.20 ±0.1 8 


0.24 + 0.13 


MS 1137.5+6625 


0.782 




7.5 ± 1.0 


8.5 ± 1.4 


5.6 ± 0.7 


0.43 ±0.1 8 


+ 0.26 


0.38 ±0.21 


0.24 ±0.13 


RX J 17 16.9+6708 


0.813 




4 6+11 

+■"-0.5 


4.3 ±0.8 


N/A 


0.37 ±0.33 


0.23 + 0.21 


0.88 ± 0.27 


N/A 


CLJ1226.9+3332 


0.890 


12.7 ±0.8 


11.5 + 0.8 


10.9 ±0.7 


9.3 ± 0.5 


0.34 ±0.11 


0.29 ±0.11 




06+ 013 



Table 7. The eight subsets of the sample used to calculate average temperature and iron abundance profiles. 



Classification of Subsets 
Subset # Redshift Temperature Cooling Number 



1 


z<0.4 


< 6.8 keV 


Yes 


6/6 


2 


z<0.4 


> 6.8 keV 


Yes 


4/4 


3 


z<0.4 


< 6.8 keV 


No 


3/3 


4 


z<0.4 


> 6.8 keV 


No 


4/4 


5 


z>0.4 


< 6.8 keV 


Yes 


2/2 


6 


z>0.4 


> 6.8 keV 


Yes 


2/2 


7 


z>0.4 


< 6.8 keV 


No 


5/3 


8 


z>0.4 


> 6.8 keV 


No 


9/7 



The two values in the number column represent the number of observations and the number of unique 
clusters in that sample, respectively. 
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Table 8. General information for each cluster: the overall (defined here to be from - 2r c ) temperature and iron 
abundance, the overall luminosity, the total number of counts in the cluster, the calculations relevant to the cooling 
time, and finally the subset as defined in Table [7] 



Temperatures and Iron abundances for each full cluster Cooling calculations for each central region Analysis Subset 



Cluster Name 


z 


Temperature 


Iron Abundance 


# of Counts 


£0-2 


Xo-0.5 


Proton Density 


Cooling Time 


Age 


Subset # 


Abell 1413 T 


0.142 


7.0 ± 0.3 


0.46 ± 0.08 


1 1 164 


5.11 


1.95 


0.019 


3.79 


8.55 


2 


Ahpll ?904 t 


0.152 


4.3 ±0.1 


86 + OS 


16382 


5. 1 7 


1 09 


0.22 


0.24 


8.44 


] 


Ahell 665 f 


0.182 


8.0 ± 0.4 


40 + 08 


12526 


3.89 


0.87 


0.012 


6 91 


8 1 1 


2 


ry Tfvno n+os^nt 

IVYV JUtJV.UTUJiU 


208 


3 6 ± 2 


If. -1-0 17 


2018 


1 96 


75 


08 


57 


7 83 


j 


Abell 773 


0.217 


8.0 ± 0.7 


0.32 ± 0.13 


8351 


9.24 


5.24 


0.0058 


14.72 


7.74 


4 


Abell 1763 


0.223 


7.9 ± 0.2 


0.34 ± 0.04 


48845 


8.81 


3.52 


0.0062 


12.93 


7.68 


4 


rtoen zjyu 


228 


8 5 ± 2 


U.4J ± U.U4 


36705 


10 70 


3 41 


038 


1 94 


7 £i 

/ ,OJ 


2 




n o<;i 

U.ZJ J 


r\ j. n i 
o.u ± U. 1 


a io ±nni 
u.jy ± u.ui 


0/10 J 1 A 
Z4ZZ 14 


1 « OA 
10. zu 


7 OA 


n Aoy 
u.uoo 


A 7 l 

u. / 1 


7 

/.Jo 


i 
i 


AKoll £Q7 
ADeil OV / 


n oqo 
u.zoz 


9.7 ± 0.7 


A 1^ _i_ a no 
U.JJ ± U.UO 


1 1/1/10 
1 J44Z 


i a a 
1 u . u 


7.75 


A AAQO 

u.uuoz 


11 i <i 
1 1 . i j 


7 1 A 
/. 1 u 


A 

4 


ADeil Oil 


a owe 

U.ZOO 


n r\ j. n i 

IX) ± U.J 


fill j.f| (\f. 

U.J 1 ± u.uo 


1 O 1 1 
1 JZI J 


U.J J 


£A 
Z.04 


n n i o 
u.uiy 


1 QO 

j.yz 


7 A/1 
/.U4 
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The column Xo-0.5 is the luminosity from - 0.5r c , and X0-2 i s the luminosity from - 2r c . The luminosities are in erg — s ! , the times are in Gyr, and the densities are in units of 
cm" 3 . Galaxy clusters with a dagger(t) next to them have confirmed radio sources likely associated with the cluster. The temperature and iron abundance have the same units as in 
Table[6] 



Table 9. Average temperature as a function of counts radius, separated by redshift, temperature and cooling times. 
Unweighted averages are in parentheses. 



Temperatures in keV 


Subset # 


Region #1 


Region #2 Region #3 


Region #4 



1 4.67(4.18) ±0.04(0.27) 6.51(6.38) ± 0.07(0.63) 7.05(6.92) ± 0.08(0.87) 7.26(7.32) ± 0.09(1.23) 

2 7.06(6.91) ±0.16(0.37) 8.37(8.04) ± 0.24(0.50) 8.68(8.36) ± 0.29(0.60) 7.51(7.45) ± 0.24(0.92) 

3 4.81(5.15) ±0.23(0.57) 5.54(5.55) ± 0.31(0.11) 5.0(5.69) ± 0.29(0.89) 4.18(5.52) ± 0.22(0.91) 

4 7.81(8.13) ±0.28(0.49) 8.69(9.12) ± 0.35(0.28) 7.34(8.61) ± 0.24(0.59) 6.90(8.15) ± 0.24(0.57) 

5 4.11(4.38) ±0.11(0.30) 5.81(6.26) ± 0.18(0.49) 7.16(7.29) ± 0.28(0.15) 6.08(5.89) ± 0.38(1.06) 

6 5.65(6.13) ±0.52(0.85) 7.68(7.95) ± 0.96(0.79) 8.12(8.27) ± 0.97(0.68) 5.56(5.55) ± 0.51(0.31) 

7 5.34(6.86) ±0.32(1.81) 4.93(4.78) ± 0.33(0.23) 4.94(4.97) ± 0.36(0.38) 4.21(2.75) ± 0.37(1.44) 

8 9.47(9.89) ± 0.37(1.09) 9.48(9.78) ± 0.37(0.55) 9.83(10.03) ± 0.40(0.63) 7.90(8.18) ± 0.29(0.63) 



Table 10. Average iron abundance as a function of counts radius, separated by redshift, temperature and cooling times. 
Unweighted averages are in parentheses. 



Iron Abundance in Solar Units 


Subset # 


Region #1 


Region #2 Region #3 


Region #4 



1 0.48(0.55) ± 0.02(0.09) 0.37(0.45) ± 0.02(0.12) 0.27(0.30) ± 0.02(0.05) 0.25(0.32) ± 0.02(0.09) 

2 0.42(0.42) ± 0.04(0.03) 0.38(0.37) ± 0.05(0.04) 0.33(0.37) ± 0.05(0.07) 0.27(0.39) ± 0.04(0.12) 

3 0.46(0.45) ± 0.08(0.07) 0.34(0.34) ± 0.07(0.04) 0.28(0.27) ± 0.06(0.05) 0.13(0.24) ± 0.06(0.10) 

4 0.44(0.43) ± 0.06(0.04) 0.20(0.25) ± 0.05(0.07) 0.37(0.47) ± 0.05(0.09) 0.27(0.25) + 0.04(0.06) 

5 0.56(0.43) ± 0.07(0.24) 0.51(0.92) ± 0.07(0.44) 0.32(0.22) ± 0.05(0.13) 0.43(0.35) ± 0.07(0.13) 

6 0.66(0.71) ±0.19(0.26) 0.26(0.26) ± 0.15(0.08) 0.28(0.31) ± 0.14(0.11) 0.39(0.45) ± 0.10(0.21) 

7 0.50(0.41) ± 0.09(0.07) 0.34(0.32) ± 0.1 1(0.09) 0.34(0.47) ± 0.10(0.21) 0.08(0.20) ± 0.15(0.20) 

8 0.47(0.47) ± 0.06(0.05) 0.25(0.24) ± 0.05(0.05) 0. 16(0.22) ± 0.05(0.05) 0.22(0.23) ± 0.05(0.08) 



